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ABSTRACT 

The observed similarities between the mass function of prestellar cores (CMF) and the stellar initial 
mass function (IMF) have led to the suggestion that the IMF is already largely determined in the gas 
phase. However, theoretical arguments show that the CMF may differ significantly from the IMF. In 
this Letter, we study the relation between the CMF and the IMF, as predicted by the IMF model 
of Padoan and Nordlund. We show that 1) the observed mass of prestellar cores is on average a few 
times smaller than that of the stellar systems they generate; 2) the CMF rises monotonically with 
decreasing mass, with a noticeable change in slope at approximately 3-5 M©, depending on mean 
density; 3) the selection of cores with masses larger than half their Bonnor-Ebert mass yields a CMF 
approximately consistent with the system IMF, rescaled in mass by the same factor as our model IMF, 
and therefore suitable to estimate the local efficiency of star formation, and to study the dependence 
of the IMF peak on cloud properties; 4) only one in five pre-brown-dwarf core candidates is a true 
progenitor to a brown dwarf. 

Subject headings: ISM: kinematics and dynamics — (MHD) — stars: formation — turbulence 



1. INTRODUCTION 

Molecular clouds (MCs) undergo a highly non-linear 
fragmentation process, even prior to the emergence of 
young stars, or in regions with no apparent star for- 
mation. The fundamental reason for their fragmenta- 
tion is the presence of supersonic turbulence that orig- 
inates at large scales from various sources, such as su- 
pernovae, spiral-arm shocks, or magneto-rotational in- 
stability. The formation of dense filaments and cores 
is the natural evolution of the intersection of ran- 
domly driven shocks in the turbulent flow. Even with- 
out self-gravity, turbulence simulations with the same 
rms Mach numbers as MCs generate density contrasts 
of many orders of magnitude, with characteristic size 
and dens ity of both filaments an d cores as observed 
m MCs (|A rzoumanian et all 120111) . Because all stars 
are born in MCs, and specifically from dense cores 
found predominantly within the densest fila ments (e.g. 
lAndre et all 12010b lArzoumanian et al.l 120111) . the frag- 
mentation by the turbulence must control the early 
phases of the star formation process. It may also con- 
trol its global properties, su c h as the star formation rate 
(|Krumholz fe McKed l200l IPadoan fc Nordlund! l20ll 
and the mass distribution of stars JPadoan fe Nordlund! 
I2002t IHennebelle fc Chabrierll2008l) . but other processes 
(e.g. stellar jets and outflows, disk fragmentation, com- 
petitive accretion, ambipolar drift) can also affect the 
mass and formation rate of stars (e.g. lAdams fc Fatuzzol 
[1991 iBonnell et al.l[200ll) . 

Most (sub-)mm or extinction surveys of star- 
forming regions have resulted in prestellar core mass 
functions (CMFs) interpreted to be consistent with 
the s tellar IMF (iMqtte et all 119981: iTesti fc Sargentl 
[19981: Uohnstone et al.l 120001 120011; iMotte et al.l 120011: 



Onishi et al.1 l2002t I Johnstone et al.l 120061 : 



2006: lAndre et al.| 120071 : lEnocli et aLl 120071 .1200 



Stankeetal 



Bell oche et al few studies have also claimed the 



detection of a CMF peak, which has been interpreted as 
a scaled progenitor of the stellar IMF peak, with the scal- 
ing mass factor giving the local efficiency of star fo rma- 
tion (lAlves et al 120071 : iNutter fc Ward-Thompson1l2007l: 
Rathborne et al.ll2009t lAndre et al.ll2010l: iKonvves et al.l 



20101 ). However, the peak of the CMF is quite close to 
the estimated completeness limits of the surveys, which 
depends also on the core shape and selection method (e.g. 
iKainulainen et al J l2009t I Pineda et al.ll2009D . 

Simple theoretical considerations show that a one-to- 
one relation between the mass of observed prestellar cores 
and that of stars is unlikely. Prestellar cores must grow 
for some time before they can collapse and form stars. 
When a core is observed it is unlikely to have just reached 
its final mass. Furthermore, the largest observable mass 
of a core may be even lower than the full prestellar mass 
involved. Once a core mass grows beyond it s Bonnor- 
Ebert (BE) mass (jBonnorl Il956l : lEbertl I1957T ) . the core 
rapidly collapses within about a free-fall time, but the 
accretion flow that was assembling that core is likely to 
still be active, and to continue to bring additional mass to 
the new protostar. The prestellar core as such is gone, 
but the prestellar mass feeding the new protostar may 
keep coming. Note that this process is essentially iner- 
tia!, and therefore distinctly different from competitive 
accretion, which assumes the flow is caused by gravita- 
ti on forces. 

iClark et"aLl (|2007[ ) have stressed a timescale problem 
in relating the CMF to the stellar IMF: if cores of dif- 
ferent masses are assumed to all contain one BE mass, 
small cores must be denser and free-fall more rapidly 
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than larger ones. If the CMF is stationary in time, 
the corresponding stellar IMF should then be steeper 
than the CMF. However, the observations do not show a 
strong correlation between core mass and density, so this 
timescale issue may not be the main problem in deriv- 
ing the IMF from the CMF. In this Letter, we study the 
relation between the CMF and the stellar I MF based on 
our IMF model (jPadoan fe Nordlundl[200l . 

2. RANDOM REALIZATIONS OF THE IMF MODEL 

The IMF model bv lPadoan fc Nordlund! (pOOl (PN02 
hereafter) predicts the mass distribution of gravitation- 
ally unstable cores generated by a turbulent flow. The 
mass of a core is the total mass, m accr , the turbulent 
flow would assemble locally, irrespective of when the 
core should collapse and cease to appear as prestellar. 
In PN02, the mass distribution of all cores (unstable or 
not) is a power law with slope x — 3/(4 — 6), where b is 
the slope of the velocity power spectrum of the turbu- 
lent flow, E{k) ~ k~ b . The Salpeter slope, x = 1.35, is 
thus recovered if b = 1.78, consistent with power spectra 
derived from the largest simulations of supersonic and 
super- Alfvenic turbulenc e and from observations (e.g. 
iPadoan et afll2006l [2001 ) . Unstable cores are then se- 
lected as those more massive than their BE mass, as- 
suming all cores have the same temperature, while their 
external density (the local postshock density) follows the 
Log-Normal pdf of supersonic turbulent flows. Most of 
the massive cores contain more than one BE mass, so the 
IMF above approximately one solar mass is predicted to 
be a power law with slope close to Salpeter's. Smaller 
cores are usually less massive than their BE mass, so the 
IMF is expected to peak at a fraction of a solar mass, 
and to decline towards smaller masses. 

This assumes that the actual stellar mass is m* = 
em accr , and the local efficiency, e, is approximately in- 
dependent of mass. In observational studies, the local 
efficiency, e corc , is instead defined as the ratio of the re- 
sulting stellar mass, m*, and the current mass of prestel- 
lar cores, m, so e cor c = m±/m. It follows that the ratio 
of the total core mass predicted by the model and the 
current core mass can be expressed as the ratio between 
our theoretical local efficiency and the observational one, 
"Wr/m = ecore/e, and thus e corc /e > 1. 

To create a population of observable prestellar cores 
consistent with the IMF model, we first generate random 
values of total core masses, m accr , following the predicted 
power law distribution with slope x = 1.35. We then 
generate a random value of external density for each core, 
according to the Log-Normal gas density distribution of 
the turbulent flow, converted to mass fraction, 



p(p)dp oc exp 



(\np + a 2 /2f 
2a 2 



dp, 



(1) 



where p is the core external density in units of the cloud 
mean density, p — p/po, and the standard deviation, a, 
of the logarithmic density field is 
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where Ms,o is the rms sonic Mach number, A^s,o = 
Vq/cs,o, with vq the three-dimensional rms velocity and 



cs,o the isothermal sound speed corresponding to the 
mean temperature To, and j3 is a characteristic ratio 
of gas to magnetic pressu re in the postshock gas (see 
IPadoan fc Nordlundl l27)lTl eqs. 16, 17, 18, 27 and 28). 
The corresponding standard deviation of the linear den- 
sity field, Up, is given by 



(3) 



Finally, we associate a random age to each core, assum- 
ing for simplicity that the star formation rate (SFR) is 
uniform over time and independent of core mass. 

In order to relate the PN02 model to observations of 
prestellar cores (cores without a detectable embedded 
protostar), we also need to account for the time evolution 
of cores, and define when, during its growth, a core would 
be observed as prestellar. In PN02, cores are assumed to 
be chunks of dense postshock filaments (or sheets), with 
size equal to the postshock thickness, A, and mass 



(4/3)7rp(A/2) 3 , 



(4) 



where p is the postshock density. Assuming that a com- 
pression from the scale £ has a shock velocity that follows 
the second order velocity structure function of the tur- 
bulent flow, v ~ £ a (with a = (6— l)/2), the compression 
lasts for an accretion time given by 



t^ CI =£/v = t (£/L ) 1 



(5) 



where to is the cloud crossing time, defined as the ratio 
of cloud size, Lq, and rms velocity, vq, to — Lq/vq. As- 
suming that the postshock thickness, A, grows linearly 
with time, we can then model the mass evolution of a 
prestellar core, up to the time t — t accr , when the total 
core mass, m accr , is reached, with the simple law 



m(t)/m a 



(^/^accr) : 



(6) 



which shows that cores must spend a significant fraction 
of their lifetime at a mass significantly lower than their 
final mass, and therefore the relation between the masses 
of prestellar cores and stars cannot be trivial. 

In order to express t acC r and m as a function of m accr , 
instead of I, we can relate m accr and t using equa- 
tion Q with the shock jump conditions, p = PqA4a,i 
and A = £/A4a,£, assuming that the shock Alfvenic Mach 
number scales like the shock velocity, M.a,i ~ v ~ £ a , 
and with t he effective rms Alfvenic M ach number ex- 
pression of IPadoan fc Nordlundl (|2011D (equations (13) 
and (14)), Ma,o = cr~. We then obtain the following 
expression for t accr expressed as a function of m accr 



mo 



(7) 



We assume that cores that do not reach their BE mass 
are seen only during their formation time, i aC cr, while 
those growing past their BE mass continue to be observed 
as prestellar for one free-fall time, reaching a maximum 
prestellar mass, m max , given by 



^be + tg 



(8) 
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Fig. 1. — Ratio of the total core mass and the current core 
mass, m. Dots are for cores that will never collapse into stars, 
plus symbols for cores that form stars and can reach m aC cr while 
prestellar, mf = m acC r < m ma x, and diamond symbols for those 
that form stars but do not reach the mass m accr while prestellar, 
m f = m max < m accr . 

where £be is the time when cores reach their BE mass, 

/ \-l/3 

*BE = *accr i (9) 

V m BE J 

and the BE mass (|Bonnorlll95l lEbertl 119571 ) is 

m BE = 1-182 cl/iG^p 1 ' 2 ), (10) 

where cs is the isothermal sound speed in the cores, cor- 
responding to the mean core temperature T, p the post- 
shock density (assumed to be the external density of the 
BE sphere), G the gravitational constant, and tg is the 
free-fall time, tg = (3ir/(32Gp)) 1 ^ 2 . The value of m max 
may be larger than m accr (if £be + tg > t accI ), and there- 
fore the final (observable) mass of a prestellar core as 
such is 

mf = min [m accr , m max ] . (11) 
3. RESULTS 

We show results for a model with characteristic molec- 
ular cloud (MC) parameters, T = 10 K, T = 7 K 
(the core mean temperature), L = 10 pc, p = 2 x 
10" 21 g/cm 3 , M s ,o = 25, = 0.4. The total MC mass 
is then to = 1.1 x 10 5 M Q , and its three-dimensional 
rms velocity vq = 4.9 km/s. We consider the simple 
case of a constant SFR for a time equal to the cloud 
crossing time, to — L /v — 1.9 Myr, and study the 
prestellar core population at a time t = to, assuming 
that prestellar cores are detected above a minimum sur- 
face density of iV"det = 10 21 cm -2 , corresponding to the 
Herschel 5-cr detect i on limit due to cirrus noise in Aquila 
([Andre et al.ll2010t lK5nvves et al.H201Q[ ). We generate 
a random distribution of total core masses, m accr , with 
probability following a power law with slope x = 1.35, 
with a minimum mass of 0.01 M©. The total mass in 
the cores is 0.26too, while the total mass of those that 
collapse into stars gives a final star formation efficiency 
of SFEf=0.05. At the time we study the core popula- 
tion, t = to, the total mass in stars (defined as all the 
cores that have reached their total mass, TO acC r) is such 
that SFE(to) = 0.02, which is a reasonable value for MCs 
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FlG. 2. — Histograms of the ratio m accr /m shown in Figure [T] 
The solid line is the histogram for all cores with mass m > 0.1 Mq, 
yielding a mean value of (Waccr/Ti) = 5.0. The dashed line in 
Figure [5] is the histogram including only core masses above half 
their BE mass, m > yielding a mean value of (m accr /m) = 

2.6. 

(e.g. lEvans etalll2009h . The core formation efficiency at 
t = t p is CFE(t n )=0-01 1 also consistent with observations 
(e.g. lEnoch et al.|[2007l ). 

Figure Q] shows the ratio between the total core mass 
and its current mass at time t = to, m accr /m — e COTe /e. 
The plot shows that prestellar cores may have a mass 
significantly smaller than the mass of the stellar system 
they will form, assuming reasonable values of e, which 
means that e coro > 1, contrary to the usual assumption 
of observational studies. Observed prestellar cores with 
masses between 0.1 and 1.0 M©, for example, may be on 
their way to form stellar systems with masses 10 times 
larger. 

The relation between current core mass and total 
core mass is quantified by the histograms of the ratio 
^accr/™, shown in Figure O The solid line is the his- 
togram for all prestellar cores with mass m > 0.1 M©, 
showing a broad distribution, with a mean value of 
(m accr /TO} = 5.0. Cores that will never grow above their 
BE mass to form stars (m accr < tube) are shown as 
dots in Figure [TJ These cores are not included in the 
histogram, but they certainly contaminate observational 
samples. 

The dashed line in Figure is the histogram including 
only core masses above half their BE mass, m > meE/2. 
These cores are somewhat closer to their corresponding 
total core masses, with an average ratio of (m accl /m) = 
2.6. The ratio of the peak of the stellar IMF from our 
model (see Fi gure [4] below) an d the peak of the multiple 
system IMF of lChabrierl (|2005D is approximately 2.1, cor- 
responding to e fa 0.48. With this value of the theoretical 
local efficiency of star formation, the average ratio be- 
tween the actual mass of stellar systems, m*, and the cur- 
rent mass of cores, to, would be e corc = e(m accr /m) = 2.4 
for prestellar cores with m > 0.1 M©, and e cor c ps 1.2 for 
the core subsample with to > tobe/2. 

Figure [3] shows the ratio of the current core mass to 
its BE mass, at t = to- It shows that a significant frac- 
tion of the cores with to > 0.1 M© may be found to 
have a mass larger than their BE mass. It also shows 
that the selection of cores with mass m > tobe/2 allows 
partial decontamination from cores that will never col- 
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Fig. 3. — Ratio of current core mass to BE mass, plotted versus 
the current core mass. The solid line shows the ratio between the 
number of cores that will never collapse and form stellar systems 
and that of true prestellar cores, after selecting cores with m > 
m BE /2. 

lapse into stars, while still providing a large enough core 
subsample to allow a meaningful estimation of the CMF. 
This strategy is already being im plemented in the ana l- 
ysis of observational surveys (e.g. IKonvves et alj I2010). 
The solid line in Figure [3] gives the ratio between the 
number of non-prestellar cores and that of true prestel- 
lar cores (computed from five different realizations of the 
same model), as a function of core mass, for cores with 
mass m > tobe/2. Around 1 M©, only 16% of these 
cores are not prestellar, while that fraction grows to 35% 
around 0.1 M Q . 

Smaller cores that may be progenitors of brown dwarfs 
can barely exceed their BE mass (towards the end of their 
lifetime). Only one in five of the selected pre-brown- 
dwarf core candidates with mass 7tibe/2 < m < 0.15 Mq 
(assuming e » 0.48) is a true pre-brown-dwarf core. Half 
of them are growing progenitors of higher mass stars, 
and one third are non-prestellar cores. Distinguishing 
true pre-brown-dwarf cores or non-prestellar cores in ob- 
servational surveys is a difficult task that may require 
the characterization of such cores using synthetic obser- 
vations of turbulence simulations. 

The CMF at the time t — t is shown in Figure |4] 
by the black histogram (and by the continuous red line 
for the average of five different realizations of the same 
model). The CMF seems to change its slope at two mass 
values, mi « 0.5 M (approximately the peak of the 
model IMF) and mi ~ 4 M Q . For masses above rri2, the 
slope is steeper than x — 1.35, it becomes shallower than 
that for masses below 1712, and it steepens again below 
mi, but still remaining a bit shallower than x = 1.35. 
The slope for masses below mi is presumably very dif- 
ficult to derive from observations, because core samples 
are generally incomplete at such low masses, depending 
on the survey sensitivity, confusion, and method of core 
selection. 

We have verified that the values of mi and m,2 scale 
1/2 

approximately as p (mi remains close to the peak 
of the model IMF), while they are approximately the 
same in models where equation (|6]) is modified to assume 
m(t) ~ t 2 , and in models where cores are assumed to re- 
main prestellar for only a fraction / of iff after they reach 
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Fig. 4. — Black histogram: CMF at time t = to. Blue histogram: 
CMF for a prestellar core subsample with masses m(fo) > »tibe/2. 
The same CMF derived from five different realizations of the same 
model is shown by the red solid line . Green h i stogra m: Model IMF 
at the time t = 2trj. Dashed lines: [Chabricr (20Q5|) system IMFs, 
shifted in mass by a factor of 2.1. 

their BE mass. However, with decreasing values of /, the 
slope of the CMF for masses in the range mj < m < m 2 
increases, becoming approximately the same as the slope 
for m < mi when / w 0.6, and approximately equal to 
1.35 when / < 0.4. 

The model IMF seen at the time t = 2i 0) when all cores 
have reached their total mass, m accr , and have already 
collapsed into stars (with SFEf = 0.05), is shown by the 
green histogra m. This model I MF is consistent with the 
system IMF of lChabrierl (|2005D . with the masses shifted 
by a factor of 2.1, shown in Figure [4] as a dashed line 
(connected to a Salpeter IMF at 2 M Q ). This shift in 
mass corresponds to e = 0.48. 

Figure 0] shows that the observed CMF should rise 
monotonically towards smaller masses, even past the 
peak of the stellar IMF. However, a careful compari- 
son between the model and the observations requires a 
detailed consideration of sensitivity, uncertainties, and 
method of core selection of each survey, besides the use 
of specific physical parameters and age for the model, 
matching those of the observed regions. Furthermore, a 
calculation of the completeness of the surveys may bene- 
fit from a model-dependent assumption on the mass and 
density distributions of small cores (near the detection 
limit), for example to estimate the confusion arising from 
projection effects. Such studies are beyond the scope of 
this Letter and should be addressed in future works. 

As mentioned in the discussion of Figure [H prestellar 
core masses would be statistically better related to 
their corresponding stellar system masses if one were 
able to select only cores with mass m > mBE/2. This 
is illustrated in Figure 0] by the CMF of such core 
subsam ple (blue histogram), which is fit by a iChabrierl 
(2005) system IMF with masses shifted by a factor of 
2.1 (dashed line), except for the largest masses, where 
the slope is steeper than Salpeter's. This CMF peaks 
at the same mass as the model system IMF, and thus 
may be used to estimate the value of e (because in this 
case e CO ro ~ e, based on the peaks alone). Its peak could 
also be used to study possible variations of e with cloud 
properties. 
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4. CONCLUSIONS 

We have computed the observable CMF predicted by 
the PN02 IMF model, assuming c haracteris t ic MC pa- 
rameters that are shown to yield a lChabrierl (|2005f) sys- 
tem IMF with e « 0.48. Our main results are the fol- 
lowing: 1) The observed mass of prestellar cores is on 
average a few times smaller than that of the stellar sys- 
tems they generate, so e corc > 1. 2) The CMF rises 
monotonically with decreasing mass, and shows a no- 
ticeable change in slope at approximately 3-5 Mq, de- 
pending on mean density; 3) the selection of cores with 
masses m > tobe/2 yields a CMF approximately consis- 
tent with the system IMF, rescaled in mass by the same 
factor as the model IMF, and therefore suitable to es- 
timate the local efficiency of star formation, e, and its 
possible dependence on cloud properties; 4) Only one in 
five pre-brown-dwarf core candidates is a true progenitor 
to a brown dwarf. 

We have not discussed the protostellar CMF, nor the 
system IMF at time t — to (when relatively massive pro- 
tostars are still growing in mass). The definition of the 
protostellar phase faces the difficulty of accounting for 
the r ate of mass transfer bet ween the core and the proto- 
star (IMcKee & Offnerll20H . However, we have verified 
that with a loose definition of protostars as all the cores 
that have reached their final prestellar mass, m max , but 
not yet their total mass, m ac cr (core ages between teE+iff 



and £ a ccr), the protostellar CMF is significantly shallower 
than the pr estellar CMF, in agreement with the obser - 
vations (e.g. IHatchell fc F"uilenl2008t lEnoch et al.ll200l . 
and the current system IMF has a high-mass slope a bit 
steeper than Salpeter's (not observationally tested yet, 
to our knowledge), due to the fact that more massive 
cores remain longer in the protostellar phase than lower 
mass ones. 

Upcoming starles s core samples fr om the Herschel 
Gould belt survey ([Andre et al.l I201Q[ ) will allow com- 
parisons with our model predictions. In this Letter we 
have only discussed the results of a very large core sample 
and a specific set of cloud parameters, without simulat- 
ing observational uncertainties and incompleteness, and 
thus setting aside a detailed comparison with observed 
CMFs for a separate work. 
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